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CLf 1 INTRODUCTION 



ABSTRACT 

The magnetohydrodynamic equations describing an inviscid fully ionized plasma in 
the vicinity of a rotating black hole are derived from a two-component plasma theory 
within the framework of the 3+1 split of the Kerr metric. Of central interest is the 
generalized Ohm's law. In the limit of quasi-neutral plasma it contains no new terms 
as compared with special relativity. Gravitomagnetic terms appear in Ohm's law only, 
if the plasma is charged in its rest frame or the MHD-approximation is not applied. 

It is argued that a relativistic single-fluid description of a multiple-component 
plasma is possible only for cold (i.e. intrinsically non-relativistic) components. As seen 
by local stationary observers, close to the horizon the electron collision time becomes 
longer than dynamical timescales, i.e. the plasma appears to behave as particles. 
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Astrophysics of black hole systems with their magnetized 
accretion discs, magnetospheres and jets has been a fasci- 
nating field of research for some twenty years. There are a 
number of papers dealing with the general-relativistic de- 
scription of those systems. The majority of the papers is 
restricted to ideal MHD of black hole magnetospheres and 
jets (Blandford & Znajek 1977; Phinney 1983; Bekenstein 
& Eichler 1985; Camenzind 1986, 1987; Takahashi et al. 
1990; Okamoto 1992; Beskin & Par'ev 1993; Fendt 1997). 
General-relativistic resistive MHD has been discussed by 
Bekenstein & Oron (1978). Khanna & Camenzind (1994, 
1996a,b), Khanna (1997) and Kudoh & Kaburaki (1996) 
have applied resistive MHD to accretion discs of Kerr black 
holes. In both ideal and resistive MHD the authors have 
used the standard Ohm's law. 

Generalized Ohm's laws in special-relativistic context 
have been derived from invariant distribution functions by 
Ardavan (1976) and by Blackman & Field (1993). Ardavan 
shows that the simple ideal MHD condition from the stan- 
dard Ohm's law breaks down in extreme regions of a pulsar 
magnetosphere. This result suggests that also around a black 
hole there should be regions in which the standard Ohm's 
law is merely a poor approximation. 

In this paper I derive the full set of MHD equations, 
including a generalized Ohm's law, for an inviscid, fully ion- 
ized plasma from a two-component plasma theory within the 
framework of the 3+1 split of the Kerr metric. I discuss two 
major issues: (i) under which circumstances is a relativis- 



tic single-fluid description of a multiple-component plasma 
possible? (ii) Does the gravitomagnetic field play a role in 
Ohm's law? 

The equations and results of Sect. 2 are general and ap- 
ply also to an electron-positron plasma. The MHD descrip- 
tion of Sect. 3 is limited to a relativistically cold electron-ion 
plasma. 

I use a metric with signs ( h + +) and set c = 1. 

Spacetime vectors and tensors are written upright. In 3- 
dimensional space vectors and tensors are denoted by ar- 
rows (e.g. v and T) and covariant derivatives are indicated 



by vertical bars, e.g. V-,/3 1 



/JV. Superscripts i, e, 



indi- 



cate the rest frames K 1 , K c and K , in which a quantity is 
defined. 



2 TWO-COMPONENT PLASMA THEORY IN 
THE 3+1 SPLIT OF THE KERR METRIC 

In this section I will first discuss the relativistic definition 
of a plasma as the centre-of-mass fluid of its two com- 
ponents and introduce the various rest frames involved. 
Then I will summarize the 3+1 split of the Kerr metric, in 
which the equations governing the individual components of 
the plasma will be written down and finally combined (in 
Sect. [|) to obtain the equation of motion, the local law of 
energy conservation and the generalized Ohm's law for the 
plasma as single fluid. 
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2.1 Basic concepts, definitions and equations 

2.1.1 Ions, electrons and plasma 

A fluid of charged particles is described by its stress-energy 
spacetime tensor, its number-flux 4-vector and its current 
4- vector J x . Thus for interacting ions (in this section 'ions' 
could mean positrons as well) and electrons we have for the 
hydrodynamic part of the stress-energy spacetime-tensor 



e c = (Pme+Pe«o J7c + coll 



= (p mx + Pi)W?Wi + plg a ' 3 + T; 



a/3 

x coll ' 



(1) 



where p mx = n x m x (l-|-II) is the total density of mass-energy 
(including specific internal energy II) and p x is the total pres- 
sure of the species 'x' as measured in the rest frame K x of 
the species (x=i,e). W x is the bulk 4- velocity of the species. 
The interaction between electrons and ions is described by 



*/3 

x coll ' 



(2) 



with the species' rest frame particle densities n x . 

The plasma is assumed to be a perfect fluid and is de- 
fined by the sums of the ion and electron tensors: 



(pL + p')W a W 8 + P 'g af> 

N a =N? + TV,? 

and 



T aB = Tf + , (3) 

(4) 



J a = Jf + J? . (5) 

Note that, from Eq. ([|), the plasma's 4- velocity, energy den- 
sity, momentum density and pressure are uniquely defined as 
centre-of-mass quantities of ions plus electrons. This means, 
in particular, that N a 7^ n p W a , where n p would be the 
number density of 'plasma particles', since there is no more 
freedom for such a definition. The consequences of this fact 
for the possibility of describing the plasma as a single fluid 
will be discussed below. 

Electrodynamic parts of T x and T are not given here, 
but will appear as Lorentz force and electric energy in the 
laws of momentum conservation and energy conservation be- 
low. 



2.1.2 Rest frames 

The plasma rest frame K is defined as the frame in which 
the plasma's centre-of-mass is at rest 

+ _ (Pmi + Pi)Vj + (Pme + PeKe + gicoll + coll _ ^ 

Pm + P ' 

(for the definition of momentum densities see below). The 
idea of a plasma implies that, in the plasma rest frame, 
both the ion-component (subscript 'i') and the electron- 
component (subscript 'e') have non-relativistic bulk veloc- 
ities v i and v e , i.e. the bulk Lorentz factors 7; = 1 and 
7c = 1 (this does not mean that the plasma components 
may not be relativistically hot). In this sense the rest frames 
of electrons K c , ions K' and plasma K are identical and 
quantities defined in K c and K 1 will be the same in K to 
high accuracy. 

The total densities of mass-energy as measured in K 

are 

£j = (Pmi + P\ Vi 2 )7i 2 + Eicoll ~ Pmi 



«i + e C ~ Pmi + Pn 



(7) 



where the approximations are due to v x 2 /c 2 <C 1 and the 
energy exchange in collisions being small compared to the 
rest mass of the particles. The latter is true for a plasma in 
equilibrium, i.e. behaving as one fluid. In this case collisions 
between electrons and ions are nearly elastic and also the 
energy flux S xcoll is negligible. The pressures p i ~ p\ and 
p c ~ Pc add up to the total pressure 



p =pi +p c 



(8) 



The charge densities are 

Pci — Pci = Zem p ca ~ p c cc = -en c p c = p ci + p cc . (9) 

2.1.3 The 3+1 split of the Kerr metric 

In the 3+1 split of the Kerr metric (see Thorne et al. 1986), 
spacetime is split into a family of three dimensional differen- 
tially rotating hypersurfaces of constant time with internal 
curvature. Those hypersurfaces can be mentally collapsed 
into a single 3-dimensional 'absolute space' in which time is 
globally described by the Boyer-Lindquist coordinate t. The 
line element of the Kerr metric in 3+1 notation is given by 

ds 2 = g a pdx a dx 13 

= -a 2 s dt 2 + h jk (dx 3 + f3 j dt)(dx k + (3 k dt) , (10) 

where the lapse function is identified with the gravitational 
redshift 



(-) 



(11) 



' FIDO £ 

r being the local proper time as measured by FIDOs (see 
below). The shift functions are the components of the grav- 
itomagnetic potential (5 

'd(j>\ _ 2aMr 

dt J fido S 2 



(3 r = f3" = , fir 



(12) 



and the components of the 3-metric h are 



hrr = -^, hee = p 2 , = J 2 , h jk = ; (J k) . (13) 

Note that h jk = g jk but h jk = g 3k + f3 3 f3 k /a 2 g . The metric 
functions appearing here are defined as 

A = r 2 ~2Mr + a 2 , p 2 = r 2 + a 2 cos 2 9 , 
E 2 = (r 2 + a 2 ) 2 - a 2 Asin 2 6> , w = (£//>) sin . (14) 

The parameters of the black hole are its mass M and angular 
momentum J, which define the Kerr parameter a = J/M. 



2.1.4 Quantities and equations in the FIDO frame 

Fiducial observers (FIDOs) who live in 'absolute Kerr-space' 
are stationary and locally non-rotating zero angular mo- 
mentum observers (ZAMOs). As seen in a global Boyer- 
-Lindquist coordinate system, ZAMOs are in shearing cir- 
cular motion around the Kerr black hole, described by their 
4- velocity 



U = aZ 1 {e t ~P) 



(15) 
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The quantities and equations describing physics in the 
locally flat laboratories of the FIDOs are obtained from 
spacetime vectors, tensors and tensorial equations by pro- 
jection along U and orthogonal to U with the projection 
operator 



ha/3 = gap + UcJJp 



(16) 



The relevant plasma quantities in the context of this paper 
are, for ions and electrons respectively, the bulk velocities 
Vi, v e and Lorentz factors 71, 7 e , the energy densities 

fx = (pmx + Px^x)7x + Excoll ~ (Pmx + Px#x)7x (17) 

and the momentum densities 

S*x = (p m x +£>x)7xi?x + SkcoU ~ (Pmx +Px)7x«x • (18) 

The stress-energy tensor for an interacting fluid in 'absolute 
space' is given by 



T*= (Pmx + Px)7x#x ® «x + Px h + Txcoll 



(19) 



The plasma can be treated as a single fluid, if collisions be- 
tween ions and electrons involve energies and energy fluxes, 
much smaller than the rest mass of the particles. Then the 

collision terms can be neglected in Eqs. ( [it] ) and ([l8|) . In T x 
the collision term has to be kept, however, as it describes 
the momentum flux, which couples ions and electrons. Note 

that e x , 5 X and T x result from splitting T x . 

The charge densities and current densities of the plasma 
components are, respectively, 



Pa = Pci7i = 



Pci 



7(1 — V- Vi) 



Pec = Pec 7c 



7(1 — V-Ve) 



(20) 



ji = PciVi j c — p cc Vc . (21) 

Equations (^) are derived in appendix [a|. 

In the context of a two-component theory, each plasma 
component (subscript x=i,e) has to obey the continuity 
equation for the density of rest-mass (see appendix [b]) 

d 



(3 ■ V (n x m x 7x) + V • {a s n^m^A) = 



(22) 



and the law of momentum conservation (Thorne et al. 1986; 
Thorne & Macdonald 1982) 



a. \ot 



V ) S x id = £x id g+ H -5xid 



1 -i 



V ■ (a g Txid) + Pcx-B + j x x B + fc x , 



(23) 



where 



fc x = e x coil g+ H ■ SxcoH V ■ (a g T xco ii) 

- — (^7 - ■ V) Sxcoll W - — V ■ (Og Txcoll) (24) 
Og \Ot / Qf g 

is an interaction force density due to collisions between ions 
and electrons. Note that, on physical grounds, ki = — k c . 
The subscript 'id' refers to the ideal (inviscid) part of the 
quantity. Equation ( pi| ) contains all the terms known from 
flat space plus an accelerating force along f3 on the l.h.s. plus, 
on the r.h.s., a coupling term between the gravitomagnetic 
tensor field 



(25) 



and Sx, which acts as gravitomagnetic force similar to the 
electromagnetic Lorentz force. The gravitational accelera- 
tion is related to the gravitational redshift by g — — Vlna g . 

The magnetohydrodynamic law of energy conservation 
for each individual fluid is (Thorne et al. 1986; Thorne & 
Macdonald 1982) 



id-'-*) 



-V • S x + 2g ■ S x + H:T X +jx • E .(26) 



The term H:T X = HijT^ 3 can be interpreted as the rate per 
unit volume at which the gravitomagnetic field does work 
on matter. 

Furthermore the electrodynamics of the plasma is de- 
scribed by Maxwell's equations and the law of charge con- 
servation (Thorne et al. 1986) 



(!-*•♦) 



Pcx + V ■ (Qgjx) = . 



(27) 



3 ONE-FLUID DESCRIPTION OF A COLD 
ELECTRON-ION PLASMA 

In this section the two sets of equations for ions (not 
positrons) and electrons will be combined to yield the equa- 
tions of MHD, which should no longer contain quantities of 
the individual components. It is argued that this is possi- 
ble only for a cold plasma. 'Cold' means intrinsically non- 
relativistic, i.e. Px "C Pmx c2 i Pmx ~ n x rn x (specific internal 
energy density negligible compared to rest mass density) and 
p m w niTTii + n m . This limit has to be taken of Eqs. (|l]) 
and (|), or, equivalently, in Eqs. fl), (@) - ©. 



3.1 Particle conservation 

Describing a plasma as a single fluid means that one has to 
define 'plasma particles' (subscript 'p'), which have to fulfil 
particle conservation 



Tjjj - ■ V ) (n P 7p) + V • (a g n p 7 P i7p) = , 



(28) 



or N°: a = (n p Wp)- a ~ 0. As was mentioned in Sect. 2.1.1, 
identifying v p with the centre-of-mass velocity v (or Wp 
with W a ) would be in conflict with Def. (|^). This means 
that, without any further assumptions, 'plasma particle' 
conservation can not be formulated in terms of v. Adding 
the continuity equations ( p2| ) (x=i,e) of the two components 
rather than the equations of particle conservation yields 

- 0- (nimi7i + n c m e 7 e ) 



+ V • (a g nimi7it/i + a s n e m e 'y B v B ) — 



(29) 



We want to eliminate the 71, y e , vi and v c in accordance with 
Definition (^) and derive a set of equations that contains 
only quantities of a one-fluid description. T Q 
splits into 



T°" 3 + T c ° 



e = e ; + e c 
and 

S = Si + S e . 

Definition (H) implies that 



(30) 



(31) 
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(Pm + i 



' 2\ 2 

v )7 



Pm7 



S : 



(Pm +P )f 2 V « Pm7 2 « 



(32) 
(33) 



like obviously a FIDO would define energy density and mo- 
mentum density of a plasma with velocity v and Lorentz fac- 
tor 7 anyway. In the cold plasma limit and with 7 C w 7 j ftt 7, 
the last five equations, together with Eqs. ( |l7| ) and (|l8|), 
yield 



(Pm7) + v • ( Q gPm7^) = 



(34) 



It is important to note that this continuity equation for 
the plasma, expressed in terms of centre-of-mass quantities, 
could only be derived in the cold plasma limit. This sug- 
gests (though a strict proof has not been given) that a hot 
multiple-component plasma may not be described as a single 
fluid at all. This has far reaching consequences for modelling 
hot astrophysical plasmas in jets and in hot transsonic ac- 
cretion flows onto black holes. 

Note that 7 C ~ 71 is not really a new assumption, but is 
a consequence of 7 C = 7i = 1. Furthermore, with the cold 
plasma assumption and nimi 3> n c m c , Eqs. ([^), ( ^ ) and 
Eqs. @, (§) imply 

7 2 « 7 ; 2 and 7 2 « w "ffvi . (35) 



3.2 Equation of motion 

In the cold plasma limit T x jd can be eliminated from 
Eq. (^3|), using Eq. (p2]), to yield the equation of motion 
for each plasma component: 



Pmx7x 



_1_ <9 

q s dt 



(7x«x) 



: Pmx7xff 



+Pmx7x H -v x V(a g p^) + pcx£ + jx x B + fc x (36) 

The equation of motion for an inviscid electron-ion plasma 
in one-fluid description follows by adding Eqs. (|36[) (x=i,e): 



Pm7 



J_ d_ 

a E dt 







(-yv) = p m 7 2 g 



+p m7 2 ff V(a gP ) + pc-B + j x B . (37) 

It does not come as a surprise that this equation looks just 
like Eq. (^6|). This is, however, not trivial and includes defini- 
tions and approximations that require thorough discussion, 
in particular in the relativistic context. This will be the topic 
of the rest of this section. Note that Eq. (^) could be de- 
rived straightforwardly by adding up the momentum equa- 
tions (see Sect. B.2.1^ b elow) and eliminating T=T\ + T c 
with the aid of Eq. (|34|) . This way, however, the physics be- 
hind the assumptions that have to be made (or that have 
been made already) would remain obscure. 

As has been elaborated in the last section, energy den- 
sities and momentum densities add according to Eq. (|^). 
Adding the current density 4- vectors of each plasma compo- 
nent and splitting into charge density and current density 
yields 



p c = p c i + Pec = Zemy, — en e7c , 
j = ji + j c = Zem^iVi — en c ^y c v c . 



(38) 
(39) 



The expressions for p c and j involve Eqs. (|^), (^o|) and 
(pl|). Note that the collision force densities have cancelled 
in Eq. @. 

Besides using the definitions above, there is a fundamen- 
tal assumption that has been made to arrive at the l.h.s. of 
Eq. (^). It has been assumed that the two plasma compo- 
nents are coupled sufficiently strong that their bulk acceler- 
ations 



rf( 7 x^x 



7x d 
a s dt 



+ 7x 



a a 



(Tx^x) 



(40) 



are nearly the same (standard assumption for the applica- 
bility of MHD), i.e. 



d(7i«i) d(y e v e 



(In 



dr c 



d("/iVi) 



(41) 



where r x is the proper time in the frame K*. The appli- 
cation of this assumption (see Schhiter 1950; Kippenhahn 
& Mollenhoff 1975 for the non-relativistic equivalent) in the 
derivation of Eq. (^) is elaborated in appendix (Q). In Kerr 
metric we have to go even a step further for consistency. The 
acceleration due to the gravitomagnetic force is of the same 
order as the acceleration along (3. In fact, if both gravito- 
magnetic terms of Eq. ( fjil) are added together their axisym- 
metric toroidal components cancel. If inequality dill ) were 
applied alone in the derivation of the generalized Ohm's law 
(see below and appendix [e]), such an axisymmetric toro- 
idal term would remain as an artefact dynamo term. We 
therefore also have to assume that the coupling between 
ions and electrons synchronizes the gravitomagnetic accel- 
erations, i.e. 



H ■(■yfvi)- H -(^v e ) « H •(7 2 w i ) 



(42) 



This assumption has no effect on the equation of motion, 
but is crucial in the derivation of Ohm's law. 



3.2.1 Exact law of momentum conservation for a 
two -component plasma 

As was mentioned above, Eqs. (j^j (with x=i,e) can be 
added to give the exact law of momentum conservation for 
a two-component plasma 

1 (^--0-V)S = eg+H-S-—V-{a s T) 



Q g \dt 



+ p c E + j x B 



(43) 



Here e, S may be the full expressions for a hot plasma. The 
stress-energy space-tensor of the plasma is 



T= (p m +p )7 v® v+p h'r 



Pm"1 2 V ( 



> v + p h 



(44) 



in consistence with Def. (|3j). The fact that, in the cold 
plasma limit, Eqs. ( |43[ ) and (^) can be combined to yield 
Eq. ( |37j ) , without explicit use of inequality ( pl| ) , implies that 
the collision terms have to implicitly guarantee synchro- 
nized bulk acceleration of the components for Def. (^) to 
be self-consistent in the first place. In the cold plasma limit, 
Eq. (^i|) can be used alternatively to Eq. (^). This is partic- 
ularly appropriate for numerical treatments of MHD, which 
requires the solutions of a conservative system of differential 
equations. 
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3.3 The law of energy conservation 

It is straightforward to show that the laws of energy conser- 
vation for the individual plasma components add linearly to 
give 



OLa \dt 



p-V)e = -V-S + 2g-S+ H:T +j ■ E . 



(45) 



Like Eq. (|43|) this equation is exact, but it makes sense as 
part of a complete set of equations only with the single-cold- 
fluid assumptions discussed above. 



3.4 The generalized Ohm's law for an electron-ion 
plasma 

Multiplying the equation of motion for the ions by Zm e 7 2 
and the equation of motion for the electrons by m^yf and 
subtracting the latter from the former yields 



2 d(7iwi 
Zy e m 



2 d(7 c « c 

7i "e" 



2 2 

= mjm e 7j 7 e 



dr\ dr c 

(Zn\ — n c )g+ H -(Zrnvi — n c v c ) 



Zm e % 



V(a g pj) 



mi7; 



V(a? g pc) 



+ (Zm e 7epci 
+Zm c '-fcki — mi7; 2 fc c 



mi7i 2 pcc)-B + {Zm e ^lji 



(46) 



This equation can be simplified by dropping terms of the 
order m e /mi as compared to their companion terms: With 
0{"fi) = 0(%) and the assumptions p\/p% <C m\/m e , 
Pci/pcc <C mi/rrie, the terms with p\, p cl and fcj can be ne- 
glected. The l.h.s. of the equation and the gravitomagnetic 
term are approximated according to inequalities (Ull) and 
(^), respectively. Finally the terrnf^) x B can be manip- 
ulated as elaborated in appendix (O) to arrive at 

m e (Zni7 2 — n c 7 2 ) / d(yv) 



en e 7e7i 



H i-yv) 



H (Zm - 



j x B V{a s p% 
en e 7 e 
1 



en 7 c a g 



nc)g 



en c 7 c 



(47) 



where r p is the proper time in the plasma frame K . 

In analogy to the classical two-component theory and 
consistent with the special relativistic treatment of Black- 
man & Field (1993) one can make the ansatz (motivated by 
the fact that k c is due to the average momentum transfer 
between electrons and ions) 



fc c = n e v c m e (-yiVi — ~/ e v e ) 



(48) 



The effective electron collision frequency v c (electron-ion, 
electron-electron, or electron-wave collisions as measured 
in the plasma rest frame) may be anisotropic. Classically 
anisotropy of u c is due to the reduced mobility of charged 
particles perpendicular to magnetic fields. It can be ex- 
pected that, as gravity, the gravitomagnetic force gives rise 
to drift motions and therefore v c might depend on the di- 
rection relative to the gravitomagnetic field as well. Addi- 
tionally, the gravitomagnetic force causes gyration of par- 
ticles around the gravitomagnetic field H = a^V x (3 



with a gyrofrequency as measured in the FIDO frame, 
~ \H\. If particles gyrate freely between collisions, i.e. if 
\H\"i/v c ~ (u}/a s )/(vc/'y) » 1 the collision rate should 
become anisotropic. A rough estimate of this ratio will be 
given below. Further discussion of an anisotropic collision 
rate is beyond the scope of this paper. 

Re-expressing j-%Vi and 7 e '?e in Eq. (^) with the aid of 
Eqs. (pi and (DS) yields 



fc e = u c m c (j - p c -yv)/e . (49) 

Inserting this into Eq. (^?]), one arrives at the Generalized 
Ohm's law 



J 
£>"7 C 



« £H —v x B — 1 £yc ' 

n e 7e en 7 e en e y e a g 

47T7o Pc-JV 

+ ^-PcffH 

4^(Zn i7 2 -n c7 2 ) (d^v) » ^ (ri(j) 



Wpe7e7i 2 



where the conductivity 



e 2 n e 



- (51) 

and oipc is the electron plasma frequency. The appearance of 
the convection current p c yv is compulsory to fulfil the trans- 
formation law of the current 4- vector and thereby confirms 
the correctness of the ansatz made for k c in Eq. (p8[). 

The generalized Ohm's law derived here contains no 
'current acceleration' term (see appendix |e|) . This is a con- 
sequence of inequalities ( f4l| ) , which has to be applied here in 
order to be consistent with Eq. (j37|), i.e. to stay in the MHD 
regime. The current acceleration term is, however, impor- 
tant only on time-scales in the plasma frame of the order of 
the collision time of the electrons, e.g. switch-on processes. 
There is, however, something like an 'acceleration current', 
which, after multiplication with o"7 e and for 7 C ~ 7i, can be 
written as 



Jacc 



v c dr u 



(52) 



A similar term can be found in Ardavan (1976) for the 
special-relativistic case and in Kippenhahn & Mollenhoff 
(1975) for the non-relativistic case. Both the 'acceleration 
current' and the 'gravitomagnetic current' 

J gm «£a£. (7 ^ (53) 

appear only in a charged plasma. 



3.4-1 The quasi-neutral limit 

In the limit of quasi-neutral plasma (Zni m n a ) and 7o « 
7i w 7 Eq. (|o|) reduces to 

j « aj(E + v x B) — (j x B) H — V(a g pl) (54) 

en e en e a g 

which contains all the terms, familiar from the non- 
relativistic generalized Ohm's law. The relative importance 
of the Hall-term can be estimated by writing 



en c c 



\j x B\ 



("c/7) 



J x 



B_ 

\B\ 



< 



("c/7) 



\J\ 



(55) 
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This shows that the Hall-term gives significant modification 
of the current, if the electron gyrofrequency measured in the 
FIDO frame, u cc = e\B\/ r ym c c, is of the order or higher than 
the electron collision rate. According to Spitzer (1962) in a 
non-magnetized plasma 



As" 



(56) 



where In A is the Coulomb logarithm. For polytropic (r = 
5/3) spherical accretion onto a non-rotating black hole (the 
solutions in Kerr metric will be similar) we have (e.g. 
Shapiro & Teukolski 1983) 



1 0.01 I rv 



T « 2 ■ 10 12 ( — I K . 



-3/2 



(57) 



(58) 



with r g = GM/c 2 and Mi = M/10M©. The dimensionless 
accretion rate m is measured in units of the Eddington rate 
MEdd = L-Edd/c 2 . Allowing for cooling, the temperature of 
a realistic accretion flow will rather be ~ 10 10 K. Thus we 
have 



and 



10 2 Mi ' 



\B\ 



\ mGauss 



/ m \ _1 /hiA\ 



(60) 



Obviously, at those high temperatures and low collision 
rates, the Hall-term would dominate the current even at low 
field strength. Note, however, that in such a situation cross- 
field driven plasma instabilities are likely to occur, resulting 
in an increased anomalous collision frequency. 

As was mentioned above, also gravitomagnetic drifts 
might cause anisotropic collision rates. The gravitomagne- 
tic gyrofrequency is given by 

(^^Mf 1 ^)^ 1 . (61) 



c 
2rT 



The ratio of gravitomagnetic gyrofrequency to Spitzer's col- 
lision rate is, as measured in FIDO frames, 



^4-3- 10 2 ^ 

K/7) 

This result hints to a problem: ui is the inverse of a dynamical 
timescale. By applying the MHD-approximation (Eqs. ^l] 
and 0|), which has led to neglecting the current acceler- 
ation term, we have implicitly made the assumption that 
our timescales of interest are much larger than the electron 
collision time. At very high temperatures Spitzer's collision 
rate is therefore too low to facilitate single- fluid behaviour of 
the plasma. Thus the gravitomagnetic field should introduce 
anisotropy of v c and of the particle distribution function. 
An anisotropic distribution function could trigger plasma 
instabilities with anomalous collision rates, which should 
make {uj / a g ) / (y c / r y) < 1. At the horizon the ratio of lo- 
cally measured gravitomagnetic gyro- to collision frequency 
blows up as 7/a g — > oo, which means that the MHD de- 
scription breaks down for FIDOs, i.e. the plasma appears to 
behave as test particles. 



3.5 Summary of equations 

As was argued above, a closed set of equations of MHD can 
only be formulated in the cold plasma limit with 7 C ~ 71 ~ 7. 
The complete set of equations for a electron-ion plasma are: 
The continuity equation 



- ■ V) (p m7 ) + V ■ {a gPml v) = 
with p m « rurrn. 

The local law of momentum conservation 



(63) 



a„ ( dt ^ y ^ 



eg+ H S V • (a g T) 

a g 

p c E + j x B , 



(64) 



with e ~ Pm7 2 1 S ~ pm'y 2 '" an( i T~ p m 7 2 iT ® v + p h. The 
generalized Ohm's law 



J_ 

(77 



+ 



E + 

4-7T7 



Zrii _ - j x B V(a g p°) 
-v x B - 1 v g ; 

n c en e j en e ja g 



Pc9 + 



Pc7« 



(T7 

47rp c / dj'-fv) 
topc'y \ d,T p 



H -(7 2 «) 



(65) 



Note that n c ~ Zp m jm\ — p c /e and p c has to be calculated 
from Eq. (Al). The electron pressure follows, e.g., with the 
assumption that electrons and ions have the same tempera- 
ture. Then p c c w Zp' /(Z + 1). 
The local law of energy conservation 



a g V dt 



V 6 = -V ■ S + 2g ■ S+ H:T +j ■ E . 



The law of charge conservation 
d 



(£-'•♦) 



P ■ V p c + V • (a g j) = . 



(66) 



(67) 



The set of equations is completed by Maxwell's equations 
(Thorne et al. 1986), from which p c and j can be computed 
without having to solve for « x in the two-fluid theory: 



V-B = 0, 

V x (a g E) = - 

V x (a g B) 



V ■ E = 47rp c , 
B . 



dt L l\ 



d_ 

dt 



E + 4:ira g j 



(68) 
(69) 

(70) 



where is the Lie-derivative along j3 , e.g. C^E = (/3 
V)E- (E- V)/3. 



3. 5. 1 Applicability of MHD and quasi-neutrality 

I have mentioned criteria that limit the applicability of 
MHD: (i) Relativistically hot plasma and (ii) if the colli- 
sion timescale of electrons becomes longer than dynamical 
timescales. The latter happens for low densities and high 
temperature, unless the collision rate is anomalous, or in a 
FIDO frame close to the horizon where y/v c — * 00. As a 
muon gas that becomes 'decayless', if it moves with high 
Lorentz factor through an observer's frame, a plasma be- 
comes collisionless in the frame of a FIDO close to the hori- 
zon. If case (ii) applies (e.g. if the Debye length becomes 
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very large) the assumption of quasi-neutrality should also 
break down. 



4 CONCLUSION 

The full set of MHD equations, including a generalized 
Ohm's law, for an inviscid fully ionized electron-ion plasma 
have been derived from a two-component plasma theory 
within the framework of the 3+1 split of the Kerr metric. It 
has been argued that the plasma can only be described as a 
single fluid, if its components are cold (i.e. specific internal 
energy density and pressure are negligible compared to the 
rest-mass energy density) . This places constraints on the ap- 
plicability of not only MHD but also of hydrodynamics to 
astrophysical objects such as hot relativistic jets and the hot 
regime of transsonic accretion flows onto black holes. Very 
close to the horizon FIDOs can no longer describe a plasma 
as a fluid, because the plasma particles become collisionless 
on local dynamical timescales. 

The generalized Ohm's law is found to be equivalent 
to the special-relativistic expression, if the plasma is quasi- 
neutral in its rest frame. Gravitomagnetic effects of the or- 
der of 'current acceleration' terms appear explicitly only in 
charged plasma, or if the MHD-approximation of synchro- 
nized bulk and gravitomagnetic acceleration of the plasma 
components is not applied. 

Unlike Ardavan (1976) and Blackman & Field (1993), 
who have derived their special-relativistic generalized Ohm's 
laws from invariant distribution functions, I have derived the 
generalized Ohm's law by combining the equations of mo- 
tion for electrons and ions, thereby compromising details of 
the interaction between electrons and ions. The ansatz made 
for the collision term is, however, confirmed by yielding, rel- 
ativistically correctly, the conduction current. 

In a companion paper (Khanna 1998) I will apply this 
theory to describe a gravitomagnetic battery, which operates 
in plasma surrounding a Kerr black hole. 
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APPENDIX A: TRANSFORMATION OF 
CHARGE DENSITIES 

Charge densities in the FIDO frame are given by p cx = 
— J^U/j,. The transformation of charge densities between 
FIDO frame and some other frame follows from express- 
ing the FIDO 4-velocity in terms of the other frame's 
4- velocity W ofJ , and its velocity as measured by FIDOs v ol i 
(see Khanna & Camenzind 1996a). The transformation be- 
tween FIDO frame and plasma rest frame K is given by 



Pcx 



7 



= ^+Jx M ^ = — +jx-«. (Al) 

7 ' 7 



The last equality follows by applying the definition of the 
current density measured by FIDOs j£ — h^J^: 



(A2) 



Here the subscript 'x' may be 'i' for ions, 'e' for elect rons 
or blank for the plasma. Inserting E ns. (|2l|) into Eq. ( Al) 
yields the second equalities in Eqs. (|20|), The transforma- 
tions between FIDO frame and electron rest frame K c , or 
ion rest frame K l , respectively, follow by inserting i7 M = 



x M /7x — v x 



into p cx 



-JJtUfi, where now x=i,e. This 



yields p cx = p? x 7x- 
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APPENDIX B: PARTICLE CONSERVATION IN 
THE FIDO FRAME 

Particle conservation is expressed in the spacetime viewpoint 
by N". a — 0, where the particle current N" = n x W^ and 
is the 4-velocity of the species 'x'. N" splits into 

n-x = h a v NZ = ?ix7x«i and Pnx = -UvNZ = n x 7 x • (Bl) 
Therefore 

NS; a = (pnxtT + <) ;ct = , (B2) 

which leads, under observance of the relation for the deter- 
minants of 4-metric and 3-metric g = — a^h and multiplica- 
tion with m x , to Eq. (p2|). 
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APPENDIX C: SYNCHRONIZED BULK 
ACCELERATION OF IONS AND ELECTRONS 



Adding Eqs. ( |36| ) (x=i,e) and using definition ( |40[ ) yields on 
the l.h.s. 

i d(7itTj) e d(j c v c ) _ ; d(jm) 

rmi j ~T~ rmc 7 /'mi 7 

ar; dr c ari 

rf(7i«i) ( d(7i«i) d(~f e V e ) 



dri y dri 



/ i , e ^ rf(Ti^i) _ ' d iiv) 

(Pmi + Pmc) j_ ~ Pn 



f/Tj 



(CI) 



which is the l.h.s. of Eq. (^). The following assumptions 
have been made in order to arrive at the last line: (i) In- 
equation ( [i~t| ) has been used and (ii) approximation ( |35| ) 
has been applied. 



APPENDIX D: MANIPULATING 

{Zm^lji - mi7; 2 J c ) 

Equation (Uq) contains the term (Zm e 7 2 ji — mHijc) x B- 
In Eq. (pO^ this term has mutated into Zenirmy 2 (jiV — 
j/Zerii) x B. This happened in the following way: The equa- 
tions 



j = ji + jc = PciVi + pceVe 

and 

-y 2 p m V W 7i 2 p mi iTi + 7epme'Ce 

can be combined to give 
Wi7i 2 J e « mi7i 2 J- Zeyrfp^v 



(Dl) 
(D2) 
(D3) 
(D4) 

PmiTi Pec 

where approximations of the order m c /mi and as given in 
Eq. (^) have been made. Putting things together and mak- 
ing the same kind of approximations yields the desired re- 
sult. 



and 



Zm e 7eji w Zp ci m e ^ 2 v 



ZuiePmelePc 



One can recognize acceleration terms, current acceleration 
terms and gravitomagnetic terms as well. If the plasma is 
quasi-neutral and for 7 C ~ 71 ~ 7, the first line and the 
second line (except j ■ V(-yv)) cancel. 

It is important to note that, any extra term containing 
j or B in Ohm's law may allow for dynamo action. In par- 
ticular, if such a term possessed an axisymmetric toroidal 
component that does not vanish in a neutral point of an ax- 
isymmetric poloidal magnetic field loop, it might support the 
growth of axisymmetric magnetic field. This would violate 
Cowling's theorem (Cowling 1934; Khanna & Camenzind 
1996a). 

The gravitomagnetic terms (/3 • V)j/en c and V/3 • 
j/en e do, in fact, individually possess such terms (±ui(jp ■ 
V lnu)e,t,/ 'eti c ], but they cancel identically. If inequali- 
ties (klh (approximation of synchronized bulk acceleration) 
were applied to Eq. (^), without approximating the gra- 
vitomagnetic force terms in the same way (Eq. |^]), the 
only gravitomagnetic term in j that remained, would be the 
V/3- j-term of Eq. (El). This would allow for an axisymmet- 
ric gravitomagnetic dynamo, but would be an artefact of an 
incomplete approximation. The generalized Ohm's law does 
therefore not introduce a new axisymmetric gravitomagnetic 
dynamo. 

The current acceleration terms are of the order 



(E2) 



An dj 

UIpe7e dr p 

They will, in general, possess axisymmetric toroidal compo- 
nents, but compared to the conduction current j/cry e they 
are of the order (^ c t p ) -1 , i.e. they violate Cowling's theorem 
only on timescales of the electron collision timescale, where 
the applicability of MHD breaks down. 



APPENDIX E: GENERALIZED OHM'S LAW 
WITHOUT SYNCHRONIZED ACCELERATION 



If the approximations of synchronized bulk and gravitomag- 
netic acceleration (Eqs. j^lj and [^2)) are not applied to 
Eq. (|l6|), and if 7 x i/ x in the bulk acceleration terms and v x 
in the gravi tom agneti c fo rce term are eliminated with the 
aid of Eqs. (D3) and (D4), one gets instead of the l.h.s. of 
Eq. ( [It] ) (or instead of the last line of Eq. [^o| with opposite 
sign) 



Zn\^ c m D 
n e eyf 



dt 



+ 



"/v -(3 



e7c 

m c 
e-y e 



To d_ 
a g dt 

1±<L + 
a K dt 



Zm 



Zm 



( 7 «) 



Zrii 



3 
eii c 



tJ2-V0. 

e7e a g 



Zrii 



(7c - 1i)v + 



(El) 
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